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Abstract 

One-dimensional Ginzburg-Landau equations with derivatives of noninteger or- 
der are considered. Using psi-series with fractional powers, the solution of the frac- 
tional Ginzburg-Landau (FGL) equation is derived. The leading-order behaviours of 
solutions about an arbitrary singularity, as well as their resonance structures, have 
been obtained. It was proved that fractional equations of order a with polynomial 
nonlinearity of order s have the noninteger power-like behavior of order a/(l — s) 
near the singularity. 

PACS numbers: 05.45.-a; 45.10.Hj 



1 Introduction 

Differential equations that contain derivatives of noninteger order ^1 |2j are called frac- 
tional equations [3J|3]. The interest to fractional equations has been growing continually 
during the last few years because of numerous applications. In a fairly short period of 
time the areas of applications of fractional calculus have become broad. For example, the 
derivatives and integrals of fractional order are used in chaotic dynamics jHlElj material 
sciences [3 |H1 E], mechanics of fractal and complex media [TUl E], quantum mechan- 
ics [121 Uni , physical kinetics [T31 [13 [HI , plasma physics [T71 [TH] , astrophysics [inj , 
long-range dissipation [201; non-Hamiltonian mechanics [211 [221; long-range interaction 
[2S1 [211 [2n] ; anomalous diffusion and transport theory [HI [201 [2Ij ■ 
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The fractional generalization of the Ginzburg-Landau equation was suggested in j2H]- 
This equation can be used to describe the dynamical processes in continuums with fractal 
dispersion and the media with fractal mass dimension [211 EOl E]- In this paper, we 
generalize the psi-series approach to the fractional differential equations. As an example, 
we consider a solution of the fractional Ginzburg-Landau (FGL) equation. We derive 
the psi-series for the one- dimensional FGL equation. The leading-order behaviours of 
solutions about an arbitrary singularity, as well as their resonance structures, have been 
obtained. 

In section 2, we recall the appearance of the Ginzburg-Landau equation with fractional 
derivatives. In section 3, the singular behavior of the FGL equation is considered. In 
section 4, we discuss the powers of series terms that have arbitrary coefficients that are 
called the resonances or Kovalevskaya exponents. In section 5, we derive the psi-series and 
recurrence relations for one-dimensional FGL equation with rational order a (1 < a < 2). 
In section 6, the example of an FGL equation with order a = 3/2 is suggested. In section 
7, the next to singular behaviour for arbitrary (rational or irrational) order is discussed. 
Finally, a short conclusion is given in section 8. 



2 Fractional Ginzburg-Landau (FGL) equation 

Let us recall the appearance of the nonlinear parabolic equation jH2l EH ESI , the 
FGL equation j2Hl 1211 EI] • Consider wave propagation in some media and present the 
wave vector k in the form 

k = ko K = ko + K|j + k;_l, (1) 

where ko is the unperturbed wave vector and subscripts (||,-L) are taken respectively to 
the direction of ko. A symmetric dispersion law uj = uj{k) for k <C /cq can be written as 



u{k) = uj{\k\) ^ u{ko) + Vg (|k| - k,) + Uj'g (|k| - ko)\ (2) 

where 



,' duj\ . ( d'^uj\ 



and 



|ko + i^\ = J(ko + + 4 ^ fco + + 7^4- (4) 
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Substitution of (jlj) into Q gives 

V v' 

Uj{k) ^UJo + VgK\\ + + if Ky, (5) 

ZKq Z 

where uq = u^ko). Expression © in the dual space ("momentum representation") corre- 
sponds to the following equation in the coordinate space: 

with respect to the field Z = Z(t,x,y, z), where x is along kg, and we use the operator 
correspondence between the dual space and usual space-time: 

u^ik) « ^ 

A generalization to the nonlinear case can be carried out similarly to (0) through a 
nonlinear dispersion law dependence on the wave amplitude: 

uj = uj{k, \Zf) ^ uj{k, 0) + b\Z\'^ = uj{\k\) + (8) 

with some constant b = duj{k, \Z\^) / d\Z\^ at \Z\ = 0. In analogy to (jHI), we obtain from 
® and O: 



(97 f)7 71 ?/ 



This equation is known as the nonlinear parabolic equation \d'2\ IS^ . The change 

of variables from {t, x, y, z) to {t, x — Vgt, y, z) gives 

that is also known as the nonlinear Schrodinger (NLS) equation. 

Wave propagation in media with fractal properties can be easily generalized by rewrit- 
ing the dispersion law @, (jH)) in the following way 



u{k, \Z\') = uj{ko, 0) + Vg^w + giii^lr^^ + (?2(4)''/' + b\Z\'\ (1< a, /5 < 2) (11 
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with new constants gi, g2- 

Using the connection between Riesz fractional derivative and its Fourier transform P] 



we obtain from (|lip 

87 87 

= -^^9^ + 9i{-^i.r'^Z + ^72(-A||)^/2Z + o;o^ + h\Z\^Z, (13) 

where Z = Z{t, x, y, z). By changing the variables from (t, x, y, z) to (t, ^,y, z), = x—Vgt, 
and using 

we obtain from (fT^ equation 
87 

I— = g,{-^^r'^Z + g2{-^\\f'^Z + cuoZ + h\Z\^Z, (15) 

that can be called the fractional nonlinear parabolic equation. For (72 = we get the 
nonstationary FGL equation (fractional NLS equation) suggested in j2H] • Let us comment 
on the physical structure of (fTH|) . The first term on the right-hand side is related to wave 
propagation in media with fractal properties. The fractional derivative can also appear 
as a result of long-range interaction |23 ■ Other terms on the right-hand-side of 
equations p3p and p5|) correspond to wave interaction due to the nonlinear properties of 
the media. Thus, equation (fT3|) can describe fractal processes of self-focusing and related 
issues. 

We may consider one- dimensional simplifications of ()15|). i.e., 

8Z 8^Z 

= + + ^^^^ 

where Z = Z(t,^), ^ = x — Vgt, or the equation 

^^ = g.^ + u.oZ + b\Z\^Z, (17) 

where Z = Z{t, z). We can reduce (fTTj) to the case of a propagating wave 

Z = Z{z-Vgt) = Z{r]). (18) 
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For the real field Z, Eq. (fTTjl becomes 

d'^Z dZ 

= ^^^^ 

where c = ivg. This equation takes the form of the fractional generalization of the 
Ginzburg-Landau equation, when = 0. 

It is well known that the nonlinear term in equations of type ()10|) leads to a steepening 
of the solution and its singularity. The steepening process may be stopped by a diffusional 
or dispersional term, i.e. by a higher derivative term. A similar phenomenon may appear 
for the fractional nonlinear equations ()19|]. 



3 Singular behavior of FGL equation 

There is an approach to the question of integrability which is not concerned with the 
display of explicit functions, but with the demonstration of a specific property. This is 
the existence of Laurent series for each of the dependent variables. The series may not 
be summable to an explicit form, but does represent an analytic function. The essential 
feature of this Laurent series is that it is an expansion about a particular type of movable 
singularity, i.e., a pole. Consequently the existence of these Laurent series is intimately 
concerned with the singularity analysis of differential equations initiated about a century 
ago by Painleve, Gambler and Garnier ISEj and continued since by many workers including 
Bureau [SZl and Cosgrove et al |HH] . 

The connection of this type of singular behavior and the solution of partial differential 
equations by the method of the inverse scattering transform was noticed by Ablowitz et al 
inn] who developed an algorithm, called the ARS algorithm, to test whether the solution 
of an ordinary differential equation was expressible in terms of a Laurent expansion. If 
this was the case, the ordinary differential equation was said to pass the Painleve test and 
was conjectured to be integrable. Under more precise conditions Conte ^H] has shown 
that the equation is integrable. Psi-series solutions of differential equations are considered 

in |iTlll2lllSllll!• 
In the paper, we consider the fractional equation 

c/D^Z(x) + cDlZ{x) + aZ{x) + hZ^{x) = 0, (20) 

where 1 < a < 2, and is the fractional Riemann-Liouville derivative: 

D:Z{x) = -J-—^ r dy—^^^^. (21) 
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Here, m is the first whole number greater than or equal to a. In our case m = 2. We 
detect possible singular behavior in the solution of a differential equation by means of the 
leading-order analysis. 

To determine the leading-order behavior, we set 

Z{x) = f{x-Xoy, (22) 

where Xq is an arbitrary constant (the location of the singularity). Then, we substitute 
fl22|) into the fractional differential equation ()20j) and look for two or more dominant 
terms. The detection of which terms are dominant is identical to the determination of 
which terms in an equation are self-similar. 

Substituting (j22I) into equation ^lU^ . and using the relation 

T(v + I) 

D>^= r( !7 (P>-1), (23) 

T[p + 1 - a) 

we get 

^/w^^^T^(^ - ^o)'"" + ^^/(^ - + - ^o)^ + hf{x - xof^ = 0. (24) 

V[p + 1 - a) 

If 1 < a < 2, then p — a < p — 1. For the dominant terms, 

^^ rdi + l-V "" " ""'^''^ ^ ^^'^"^ ~ ^ °' ^^^^ 

As a result, we obtain 

p — a = 3p, (26) 
T{p +1- a) 

Equation gives 

P=-2- (28) 
If 1 < a < 2, then —l<p< —1/2. Therefore the leading-order singular behavior is 
found: 

and the singularity is a pole of order a/2. Evidently our psi-series starts at (x — Xq)~"^'^. 
The resonance conditions and psi-series is considered in the next sections. 

As a result, we get that fractional differential equations of order a with polynomial 
nonlinearity of order s have the noninteger power behavior of order a/(l — s) near the 
singularity. 
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4 Resonance terms of FGL equation 



The powers of {x — xq) that have arbitrary coefficients are called the resonances or Ko- 
valevskaya's exponents. To find resonance, we consider the substitution 

Z{x) = f{x-Xoy + l{x-XQy+', (30) 

and find the values r. In equation (|Hn|l parameters p and / are defined by 



a / 9Til~a/2) 

^=-2' ^=V"^r(l-3a/2)- ^^'^ 

Substitution of Eq. (jSOI) into (jUl) gives 

9ff^^^^ix - ^o)"-" + cpf{x - xof-' + af{x - xof + bf{x - xof^+ 

r(p + r + 1 — a) 

+3bl^f{x - Xo)'P+^'' + 3blf{x - XoY^^'' = 0. (32) 
Using equation (jHH), and considering the linear with respect to / terms of (jH^ . we have 

r(l + r-a/2) _ , r(l-a/2) ^ 

r(l + r-3a/2) r(l - 3a/2) ^ ^ 

This relation allows us to derive the values of r. Equation ()33|) can be directly derived 
by using the recurrence relations. In the general case, the values of r can be irrational or 
complex numbers. The solution of the FGL equation with 1 < a < 2 have two arbitrary 
parameters. Therefore, we must have two values of r that are the solutions of equation 
()33|) . It is interesting to note that ()33|) gives two real values of r only for 

a > ao, (34) 

where 

ao ~ 1.3005888986. (35) 

The order ao is an universal value that does not depends on values of parameters g, a, b, 
c of the FGL equation 
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The plots of the function 

y^"^^- T{l + r-3a/2) ^r(l - 3a/2) ^^^^ 
are shown in figures 1 and 2. The solutions of equation correspond to the points of 
intersections with the horizontal axis. 

As a result, the nature of the resonances is summarized as follows: 

1) For a such that 1 < a < ao the values of r are complex or r < —a/2. 

2) For a such that ao < a < 2, we have two the real values of r. Note that for ao < a < 
1.9999999995, the values r satisfy the inequality r < 6.426. 




0,2 0,4 0,6 0,8 1 

r 

Figure 1: Plot for the order a = 1.30. 



5 Psi-Series for FGL Equation of Rational Order 

Let us consider the psi-series and recurrence relations for the one-dimensional FGL equa- 
tion (|^n|). where the order a is a rational number such that: 

Tn 

a = — , (1 < a < 2). (37) 
n 

Following a standard procedure jH], we substitute 

^ oo 

Z{x) = -^--j-,Y.e,<j>\x-xo), (38) 




Figure 2: Plot for the order a = 1.31. 



where 

into the fractional Ginzburg-Landau equation (f^n|) . Note that the coefficient Cq is a real 
number for two cases: (1) g/b > and 1 < a < 4/3; (2) g/b < and 4/3 < a < 2. 
For the rational order a = m/n, we suggest to use the 0-function in the form 

0(a;-a;o) = (x-xo)^/^". (40) 

Then 

^ OO CO 

Z{x) = _ Y.^u{x- xof^ = J2 ek{x - xof^-'\ (41) 

^•^ k=0 k=0 

where 

A = I- («) 

In this case, the action of fractional derivative of order a = m/n can be represented as 
the change of the number of term Ck ek-2m- 

D:{x - xor^ = - xo)^--. (43) 

J- {Pk-2m + J-j 
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It allows us to derive the generalized psi-series solutions of fractional Ginzburg-Landau 
equation. 

Substitution of the series 



oo 

k — m 



Z{x) = ^kix - xo)^-"/^ = ^ e,{x - xo)^ (44) 

A:=0 A:=0 

into equation gives 

oo -pf k-m+2n \ r 
fe=0 ^ 2n I fc=0 

oo CXD CX> oo 

efc(x-Xo)2n + 5 2^ 2^ 2^ efciefc^efcglx - xo) 2" =0. (45) 

fc=0 fci=0 fc2=0 fc3=0 

Let us compute (A; = 1,2,...) through the equation of coefficients of like power of 
(x — Xo) to zero in (jl^ : 

oo p/'fc-m+2nN oo , „ „ 

efc . fc_3^+2n N (^-^0) 2" +C 2^ efc_2™+2n ^^^^ (x - Xq) 2» + 

fc=0 ^ 2n k=2m-2n 

oo oo 

E fc — 3m fj ^ > fc — 3m 

efc-2m(x-xo) 2n 35eQ2_^efc(x - Xo) 2" + 

fc=2m fc=0 
oo fc— 1 j 

+6 ^ ^ ^ ek-i-jCiCj (x - Xo) = 0. (46) 

fc=0 4 = 1 j = l 



Using el = p, we get 



■^^^^^) fc-3m + 2n 



•p^ fc— jm.+2n r^'^ 2^ 

fc-1 i 

+aefe_2m + + ^ X] X] efc-i-jeiCj = 0. (47) 

i=i j=i 

Here = 0, 1, 2, and = for k < 0. We can rewrite the recurrence relation (j47|) as 

/ pj- fc-m+2n ^ „,,2l fc - 3m + 2n .v^v^ 

^fc I ^ p( fc-3m+2n -) + ^^Z 1 = -cefc_2m+2n ^^^^ aCfe.a^ " ^ efc-i-je^ej, 

(48) 
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where 

^ ^r(i - C./2) ^ _^r(^ 

6r(l -3a/2) 6r(^^)' ^ ^ 

Substitution of equation into equation PH|) gives 

r(^^f±g^) _ 3r(g^) \ _ _ A:-3m + 2n _ 

^'=5' I p^ fc-3m+2n ^ p^ 2n-3m ^ I CeA;_2r)i+2n 2^ 0,^k-2m 

k-1 i 

-h^^ek-i-jCiep (50) 
Note that we get resonances for the k that satisfies the condition 



r(^^i^) 3r(gM) 

fc— 3m.+2ra N 2n— 3m N 

2n I V 2n / 



0. (51) 



In this case, the coefficient Ck can be arbitrary. 
As a resuh, we obtain the nonresonance terms 

cek-2m+2n ' h aek-2m + ^ Ck-i-jCiej I , (52) 

''^ i=l j=l 



where 



-p^ fc— 3m+2n ^-p^ 2n— 3m 
A{k,m,n) — ~ ^^-^^ k-m+2n \-r^ 2n-3m \ f k-3m+2n \-r ^ 2n-m \] " (^'^) 



6 Fractional Ginzburg-Landau Equation with a = 1.5 

Let us consider the FGL equation with derivative of order a = 3/2. In this case, 
n = 2, m = 3, and the coefficients Cfc are defined by 



(ce,_2^ + aefc_6 + b El"/ E;=i e.-i-,e,e,) r(^)r(f ) 
For /c = 0, we have 



where r(i)/r(^) > 0, and we suppose g/b < 0. For k = 1, equation leads to ei = 0. 
For A; = 2, 



3ceor(f)r(f) 

4^7[r(|)r(f)-3r(f)r(i)]- 



^2 = ..^T..3^r^-5^' orZ-SNr^iM ' ^^6) 
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Using the relation we get 



^ cv/^5(^7r3/^2^/^r(3/4) 

2(7[2r4(3/4) + 57r2] ' ^ ^ 



For = 3, 63 = 0. For = 4, 



^ (ce,^ + feoe^)r(^)r(f) 

^7[r(|)r(f)-3r(-i)r(i)]- 



Substitution of (jHSj), and (jFfj) into (j3Hj) gives 



_ c'^5gnV2/bT\3/A) 

~ 4^2[2r4(3/4) + 57r2]2 ^^^^ 
For = 5, we have 65 = 0. For k = 6, equation is 

^ (ce4| + aeo + fej + 2feoe2e4) r(l)r(^) 

^?[r(i)r(f)-3r(i)r(i)] • ^''^ 

Substituting eo, 62 and 64 from equations (j^ . (jSTj), and into (jUUj) . we have 



+20c^r^(3/4)7r^ + 16a/r^2(3/4)^^20a/r^(3/4)7r2 + 300a/r^(3/4)7r^ + 250ac/Vj. (61) 
As the result, we obtain 

00 

= ^efc(x - xo)~ = 

k=0 

= eo{x - xo)-^/^ + 62(2; - xo)-^/^ + 64(0; - xoY^^ + eeix - Xof^^ + ... (62) 

This is a power psi-series that present the solution of FGL equation of order a = 1.5. 
The coefficients in equation (jH^ are defined by equations (jFfjl . (jSHl), and (jHT|l . For 

example, a = —b = c = g = l gives 

eo ^ 0.9615539375, 62 ^ -0.2382246293, 
64 ^ 0.001685563496, eg ^ 0.3872134448. 
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7 Next to Singular Behavior 



In sections 5 and 6, we consider the rational a. In this section, we suppose that the order 
a is an arbitrary (rational or irrational). Instead of imposing a series commencing at the 
power indicated by the singularity found by the leading-order analysis, we can determine 
the next to singular behavior by writing 

Z{x) = j{x-x^Y^G{x), (63) 

where 

We can always write Z{x) in the form ()fi3j) . To make the process useful, we require that 
the first term be the leading order term, i.e., 

(x - xo)~^G(x) = {x- xoT^'^G{x) ^0 if (x - xq) 0. (65) 

Substituting (jH!^ into (j^ . and using we have 



9f 



T{p+1 



T{p+l-a 



^-{x - xof-" + cpfix - xof-' + af{x - xoY + hf{x - xof'+ 



+gD'^G{x)+cDlG{x)+aG{x)+hG'^{x)+?,hf{x-xof^G{x)+^hf{x-xoYG'^{x) = 0. (66) 
Equations and give 

gD^G{x) + cDlG{x) + aG{x) + bG^{x) + 3bf{x - xo)^^G(x)+ 

+3bf{x - Xo)PG\x)cpf{x - XoY'^ + af{x - Xof = 0. (67) 

Multiplying this equation on (x — Xo)~^^, and using condition we get the equation 
without nonlinear terms for (x — Xq) — > 0. 
As the result, we obtain 

gD^G{x) + cDlG{x) + aG{x) + cpf{x - Xo^^^^^ + af{x - Xo)~°/^ = 0, (68) 

with condition (|65|) for the solutions. Equation (|68|) is a linear inhomogeneous frac- 
tional equation. The solution of this equation allows us to find the solution of the one- 
dimensional FGL equation. 
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Let us consider equation (jHH|l with c = 0, and the boundary conditions 

{d:-'G{x))^^^^ = G,, k=l,2. (69) 

Then the solution is 

2 

k=l 

+af [ {x-xo- yr-'E^,4-a{x - xq - yr]{y - xo)-"^^dy. (70) 
Here Ea,i3[z] is a Mitag-Leffler function that is defined by 

oo ^ 

Let us consider the integral representation for the Mittag-Leffler function 

where Ha denotes the Hankel path, a loop which starts from — oo along the lower side of 
the negative real axis, encircles the circular disc |^| < Iz]^^"" in the positive direction, and 
ends at — oo along the upper side of the negative real axis. By the replacement — > ^ 
equation (f7^ transforms into EH] : 

Ea,M = TT— / (K a < 2), (73) 

where na/2 < 5 < min{7i, na}. The contour 7(0, S) consists of two rays S^s = {arg(^) = 
—6, 1^1 > a} and S+s = {arg(^) = +6, \^\ > a}, and a circular arc Cs = = 1, — 5 < 
arc{^) < 6}. Let us denote the region on the left from 7(0, 6) as G~{a, 6). Then 



n=l ^ ' 

and 6 < |arg(2;)| < n. In our case, z = —a{x — Xq)", arg(z) = vr. As a result, we 
arrive at the asymptotic of the solution, which exhibits power-like tails for x — ^ 00. These 
power-like tails are the most important effect of the noninteger derivative in the fractional 
equations. 
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8 Conclusion 

In this paper, we generalize the psi-series approach to the fractional differential equations. 
As an example, we consider the fractional Ginzburg-Landau (FGL) equation [2H1 123 EOl 
IHT] . The suggested psi-series approach can be used for a wide class of fractional nonlinear 
equations. The leading-order behaviours of solutions about an arbitrary singularity, as 
well as their resonance structures, can be derived for fractional equations by the suggested 
generalization of psi-series. 

In the paper, we use the series 



where k, m, n are the integer numbers. For the order a = m/n, the action of fractional 
derivative 



can be represented as the change of the number of term ek-.2m in (f^ . It allows us 
to derive the psi-series for the fractional differential equation of order a = m/n. For the 
FGL equation the leading-order singular behaviour is defined by power that is equal to 
the half of derivative order. 

A remarkable property of the dynamics described by the equation with fractional space 
derivatives is that the solutions have power-like tails In this paper, we prove that 
fractional differential equations of order a with a polynomial nonlinear term of order s 
have the noninteger power-like behavior of order a/{l — s) near the singularity. 

It is interesting to find barriers to integrability for fractional differential equations. In 
general, the integrability of fractional nonlinear equations is a very interesting object for 
future researches. It has many problems that are connected with specific properties of 
the fractional calculus. For example, we must derive a generalization of the Lie algebra 
for the vector fields that are defined by fractional derivatives. For this generalization, the 
Jacobi identity cannot be satisfied, and we have a non-Lie algebra. The definition of such 
"fractional" Lie algebra is an open question and cannot be realized by a simple way. To 
formulate the fractional generalization of a Lie algebra for derivatives of noninteger oder, 
we can use the representation of fractional derivatives as infinite series of derivatives 
of integer orders fT^. For example, the Riemann-Liouville fractional derivative can be 




(75) 




(76) 
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represented as 

D: = J2M^,^o,(y)^, (77) 

n=0 

where 

A„(..a..)^ (-f-'°^'"-°)(;rrV P8) 

i (1 — ajL [n + Ij L [n + I — a) 

Then the possible reahzation of the generahzation is connected with the special algebraic 
structures for infinite jets [IHj. These structures and approaches can help to solve some 
problems that are connected with the integrability of fractional nonlinear equations. 
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